From quantum graphs to quantum random walks 
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We give a short overview over recent developments on quantum graphs and outline the connection 
■ between general quantum graphs and so-called quantum random walks. 

o\ . 

(N . 1 Introduction 



The study of quantum mechanics on graphs has become an important tool for investigating the influence 
of classical dynamics on spectra, wavefunctions and transport properties of quantum systems. Quantum 
, networks have been used with great success to model quantum phenomena observed in disordered metals 

£NJ ' and mesoscopic systems (Chalker and Coddington 1988); typical behaviour found in diffusive systems 

such as localisation - delocalisation transitions (Freche et al 1999), transport properties (Pascaud and 
Montambaux 1999) and quantum spectral statistics (Klesse and Metzler 1997) have been studied on 
graphs in the limit of infinite network size. Kottos and Smilansky (1997, 1999) looked at quantum graph 
models for general, non-diffusive graphs; this approach was motivated by trying to understand the validity 
of the Bohigas-Giannoni-Schmit (BGS) conjecture (Bohigas et al 1984) in terms of periodic orbit trace 
formula. The conjecture relates the properties of the classical dynamics of a systems to the spectral 
statistics of its quantum counterpart and states that for chaotic systems the statistics depends only on 
the symmetries of the problem and follows random matrix theory (RMT) otherwise. 

Making use of the fact that periodic orbit trace formula are exact on quantum graphs and that there 
are only a finite number of different length scales on a finite graph, a series of remarkable results have 
been obtained over the last couple of years. A closed form quantisation conditions in terms of periodic 
orbits has been given by Bliimel et al (2002); Barra and Gaspard (2000) derived an integral expression 
for the level spacing distribution starting form the periodic orbit trace formula. Furthermore, Schanz and 
Smilansky (2000) described localisation on one-dimensional chains in terms of combinatorial expression 
using periodic orbits. The most far reaching development is maybe due to Berkolaiko et al (2002,2003), 
who, inspired by work from Sieber (2002) and Sieber and Richter (2001), derived next to leading oder 
terms for the formfactor. Extensions to all orders have been given by Miiller et al . (2004). Gap conditions 
given by Tanner (2001) and Gnutzmann and Altland (2004) give lower and upper bounds for the border 
of universality. For a recent review on quantum graphs see Kuchment (2002). 

Quantum dynamics on graphs became an issue also in the context of quantum information. Aharonov 
et al (1993) pointed out that a random quantum walk on one dimensional chains can be faster than 
the corresponding classical random walk. Since then, a whole field has emerged dealing with quantum 
effects on graphs with properties superiour to the corresponding classical operations. For an introductory 
overview and further references, see Kempe (2003). 

We will in the following give a general definition of quantum graphs and discuss a specific set-up 
considered by Kottos and Smilansky (1997). We will then review recent developments on the spectral 
statistics of quantum graph ensembles. Next, we discuss a special class of quantum graph ensembles, so- 
called regular quantum graphs. These types of graphs can show strong deviations from RMT depending 
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on topological properties imposed on the graph in form of edge-colouring matrices. We will show that 
such graphs can be interpreted as realisations of quantum random walks on graphs. 

2 Quantum graphs - a brief review 
2.1 Quantum graphs on line graphs 

In its most general form, a quantum graph is defined in terms of a (finite) graph G together with a 
unitary propagator U; it describes the dynamics of " wavefunctions" <f> on the graph according to 

4>n+i = U0„ , 

such that waves can propagate only between connected vertices. Motivated by physical application we 
will adopt a construction of quantum graphs in terms of so-called line-graphs as explained below. 

A (finite) directed graph or digraph consists of a finite set of vertices and a set of ordered pairs of 
vertices called arcs. We denote by V G and E G the set of vertices and arcs of the digraph G, respectively. 
Given an ordering of the vertices, the adjacency matrix of a digraph G on n vertices, denoted by A , is 
the (0, l)-matrix where the ij-th element is defined by 

a g . = { l if(u)e£ G , (1) 

y ' |^ otherwise. 

An undirected graph (for short, graph) is a digraph whose adjacency matrix is symmetric. The undirected 
connections between vertices are called edges in this case. The line digraph of a digraph G, denoted by LG, 
is the graph which is obtained when taking the arcs as the new vertices; it is thus defined as V LG = E G 
and, given (hi),(jk) <E E G , the ordered pair ((hi)(jk)) £ E LG if and only if i = j (Bang- Jensen and 
Gutin 2001). 

A quantum graph associated with a digraph G on n vertices can then be defined in terms of a set of 
unitary vertex scattering matrices cr^) on vertices j — 1, ...n and a set of arc lengths defined for 
every arc (ij) € E G . Waves propagate freely along the directed arcs, transitions between incoming and 
outgoing waves at a given vertex j are described by the scattering matrix o"W, see Fig.^,. The two sets 
specify a unitary propagator of dimension ue = \E G \ defining transitions between arcs (ij), (i'j') € E G 
which has the form (Kottos and Smilansky 1997) 

S G = DV with D( im , r) = 6 i>if 6 jtf e ikL » , (2) 

where A; is a wave number and 

The local scattering matrices describe the underlying physical process which may be derived from 
boundary conditions on the vertices. We will construct an example below but may often regard the 
cr^'s as arbitrary unitaries. Let d~ and df be the number of incoming and outgoing arcs of a vertex 
i, respectively. A sufficient and necessary condition for a digraph G to be quantisable in the way above 
is then, that for every vertex i S V G , df — d^ = di = dimcr^ (Pakonski et al 2003). This means in 
particular that if G is an undirected graph then it is quantisable. 

Kottos and Smilansky (1997) considered solving the Id Schrodinger equation on an undirected graph 
assuming free propagation on the arcs and imposing continuity and flux conservation at the vertices. The 
solution on each arc propagating from vertex i — > j takes on the form 

<t>( Xij ) = 4,+ e ikx » = <j>T. e - ik ^-^\ 

with (f>ij being the outgoing (+) or incoming (— ) wave at vertex i or j. Continuity and flux conservation 
can then be written in terms of the amplitudes <^t-(0) and <f>~(Lij) at the vertices, that is 
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Figure 1: a) Quantum graph; b) periodic path on a graph with quantum weights defined with respect to 
the line - graph. 



Continuity: 4>% = ^ji = °i f° r au 3 with(z, j) G i? G 
Flux cons.: ^ <f>± = ^ &ji 

These conditions give rise to local scattering matrices a 1 mapping amplitudes (j)^ onto <j)^ at vertex I 
having the form 

The eigenvalue condition is then given as 

det(I-S G (/fc)) =0 (5) 

with S G (fc) = D(fc)V as defined in 10), ©• Scatt ering matrices for more general boundary conditions 
can be found in Kottos and Smilansky (1999). 

The "classical" dynamics corresponding to a quantum graph defined by a unitary propagator S is 
given by a stochastic process with transition matrix T defined by 

S G ^T(S G ) with T(S%HS G | 2 HV^-| 2 . (6) 

The matrix T is clearly stochastic, as = 1 diic to the unitary of S G ; the set of transition 

matrices related to a unitary matrix as defined in (JHJ is a subset of the set of all stochastic transition 
matrices, referred to as the set of unitary- stochastic matrices. The topology of the set in the space of all 
stochastic matrices is in fact quite complicated, see Pakohski et al (2001). In what follows, we will only 
use that T has a largest eigenvalue 1 with corresponding eigenvector ^-(1, . . . , 1) which follows from the 
Frobenius-Perron theorem and T being unitary-stochastic. 

Note that both the quantum mechanics as well as the associated stochastic dynamics relates to 
transitions between arcs and is thus defined on the line digraph of G. 

2.2 Unitary stochastic ensembles and spectral statistics 

Inspired by the BGS conjecture for quantum systems, we expect a link between the dynamical properties 
of the stochastic process T and the statistical properties of the spectrum of associated unitary matrices 
S G . It is thus natural to consider the ensemble of unitary matrices 

USE To = {S G is quantum graph on G \ T(S G ) = T } (7) 



3 



for a given graph G and a fixed unitary stochastic matrix To associated with a stochastic process on the 
line-graph of G. Clearly, if S G (fc) £ U SEt for k = then it is for all k £ 1R. In fact, if all the arc length 
Lij are incommensurate, D(fc) sweeps out the space of unitary diagonal matrices of dimension tie- For 
practical purposes we will thus often replace averages over a given U SEt by averaging over the space of 
diagonal unitary matrices D with Dij = Sije llfj using the Euclidean measure on the tie - torus. 

Rather than looking at the spectrum obtained from the secular determinant © , we will here consider 
the spectrum S G for fixed wavenumber k and than average over k. One can write the spectrum in terms of 
a periodic orbit trace formula reminiscent to the celebrate Gutzwiller trace formula being a semiclassical 
approximation of the trace of the Green function (Gutzwiller 1990). We write the density of states in 
terms of the traces of S G , that is, 

TLB OO 

d{6) = V 8(6-0i) = ^ + -Re V Tr(S G )"e- me , (8) 

* — ' 271" 7T * — ' 

i—l n— 1 

where {#i};=i,n B refers to the eigenphases of S G . The traces Tr(S G )™ can be given as sum over all 
periodic paths of length n on the graph, i.e. 

(») 

Tr(S G )™ =y A p e lkL ". 
p 

Describing a given periodic path in terms of its vertex code (vi, vi . . . v n ), V% € {1, 2, . . . n} with (vi, Uj+i) € 
E being an allowed transitions between vertices, one obtains for the amplitudes A p and lengths L p 

n n 

= JJ V(viVi + i),(v i+ iVi = 2)i L p = y^ L ViVi+1 . (9) 

i=l i=l 

where the amplitudes give again transitions between arcs (not vertices) . An example of such a periodic 
path is given in Fig. ^d. 

Trace formulas like (JSJ are a starting point for analysing the statistical properties of quantum spectra. 
The statistical quantities such as the two-point correlation function can be written in terms of the density 
of states d(8,N), that is, 

R2{x) = =f < d(0)d(0 + x/d) > US E T ,e , (10) 
d 

where d = N/2n is the mean level density. The average is taken here over the angle 9 as well as over the 
USE (which is equivalent to energy averaging). The Fourier coefficients of (1101) can be written in terms 
of the traces of S; one obtains 

K(t) =< ^|TrS^| 2 >use t (11) 

with t = n/N and the average is taken over a USE. The so-called form factor K(t) can thus be written 
as a double sum over periodic paths on the graph 



, ( n ) 

K ^ = <^J2 A P A P' eiiLp ~ Lp ' ]> usE T (12) 
p,p' 

(n) 

w <?^TrT"+ < y ApA p ,e^- L ^ > USEt ■ (13) 
p¥=p' 

The first term in (|13|l . also called the diagonal term (Berry 1985), originates from periodic orbit pairs 
ip,p') related through cyclic permutations of the vertex symbol code. There are typically n orbits of 
that kind and all these orbits have the same amplitude A and phase L. The corresponding periodic orbit 
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Figure 2: Periodic orbit with one intersection on undirected graphs: two different paths exist having the 
same length following either the black or blue orbit. 

pair contributions is (in general) g ■ n - times degenerate where n is the length of the orbit and g is a 
symmetry factor {g = 2 for time reversal symmetry). 

Expanding the random matrix result for the formfactor, one obtains for < r = n/N < 1 

Kcue(t) = t 

K C oe{t) = 2r-rlog(l-2r) 

= 2t-2t 2 + 2t 3 + ... . 



The linear terms are reproduced by the diagonal contribution which gives the important link between 
the stochastic dynamics and the spectral statistics; it follows from lim„_» 00 TrT" = 1 given there is a gap 
between the leading and next-leading eigenvalue of T. Contributions to the double sum in which 
survive the ensemble average are due to periodic orbit pairs where orbits visit the same arcs along its 
path but in different order. The simplest examples are "figure eight" orbits of the type shown in Fig. [21 
which arise in undirected graphs. Berkolaiko et al (2002) could indeed show that orbits of that type give 
the correct 0(t 2 ) - contributions to the GOE - form factor. In fact, a general scheme emerges relating 
0(r n ) contributions to periodic orbits with n intersections. This work was inspired by a similar analysis 
for general quantum systems by Sieber and Richter (2001). The periodic orbit contributions giving the 
C(r 3 ) have been worked out by Berkolaiko et al (2003) and a general scheme for obtaining higher order 
terms iteratively has been given by Muller et al (2004). 



2.3 Border of universality 

In the light of these recent developments, it becomes important to establish the boundaries at which 
the spectral statistics of quantum graph ensembles (or more general quantum systems) starts to deviate 
from random matrix behaviour. One can distinguish two different scenarios. Firstly, the properties of 
the underlying dynamics, that is, the stochastic process T may be linked to the spectral statistics and 
may thus provide conditions for the onset of deviations from RMT behaviour; such an approach is in the 
spirit of the original BGS - conjecture making a connection between classical chaos and random matrix 
statistics. Secondly, one may consider special phase or length correlations in the quantum graph which 
could lead to interesting non-universal statistics; in this approach, quantum graph ensemble averages are 
carried out only over subsets of the full USE. 

We will discuss known bounds on the border of universality related to the properties of the stochastic 
process. An interesting family of quantum graphs which belongs to the second category are so-called 
regular quantum graphs, which will be treated in more detail in section |21 



Spectral gap conditions When studying the border of universality, we always need to consider the 
limit of large graphs, that is, ue — > oo. This limit is in general not well defined, but may often be 
obvious from the examples considered. We will thus define the semiclassical limit loosely via a family 
of unitary-stochastic transition matrices {T„ E } and associate USE's and take tie — * oo. The leading 
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Figure 3: Form factor for a) star graphs: numerics (red) versus the power series expansion <|15[) (green); 
b) diffusive network (in 2 d): deviations occur for small r, see inset. 

term in l|13|) then gives a condition for a family to show deviations from RMT statistics in terms of the 
spectrum of T; the diagonal term must obey 

1 - TrT ,lBT w e - AnET -> 0; for n E -> oo 

in order to match the leading coefficient in the expansion of the RMT form - factor; here A is the spectral 
gap, that is, A = — log(l — |Ai|) with Ai being the next to leading eigenvalue of T nj5 ) and r = n/riE 
is fixed. The condition above implies that we expect to see deviations from RMT behaviour to leading 
order if 

A ~ n E a with a > 1, 

that is, whenever the gap closes faster than 1/ue for large system sizes (Tanner 2001). Based on super- 
symmetric techniques, Gnutzmann and Altland (2004) could give a lower bound by showing that the 
spectral gap condition guarantees RMT behaviour for a < 1/2. The border of universality must therefore 
lie in the range 1/2 < a < 1. Other bounds have been given by Berkolaiko et al (2002,2003) which have 
been derived from higher order terms in the expansion of the form factor. Below we we will give two 
examples families with a = 1 displaying critical behaviour by neither converging to RMT nor to Poisson 
statistics in the semiclassical limit. 

Quantum star-graphs: Quantum star-graphs arise naturally when one quantises a graph with a single 
central vertex attached to tie undirected edges, see inset of Fig.[3]j,. The underlying graph is complete, 
that is, we can reach every edge from every other edge through the central vertex (ignoring the trivial 
dynamics on the outer vertices). Typical boundary conditions imposed on the wave equation at the 
central vertex following Kottos and Smilansky's approach in Sec. 12.11 yield scattering matrices which 
greatly favour backscattering. The vertex scattering matrix Q which is essentially equivalent to the 
matrix S here is 

S ti = -6a + . (14) 

TIE 

The Markov processes associated with quantum star graphs correspond to systems of weakly coupled 
edges. Its dynamical properties are determined by the spectrum of the stochastic matrix associated with 
<|14|) which is highly degenerate and can be given explicitly (Kottos and Smilansky 1999), that is, 

4 4 

A = l, Ai,...,A nE _i = (1 )» — . 

TIE TIE 

Quantum star-graphs have therefore a critical classical spectrum with a spectral gap vanishing propor- 
tional to 1/ue', one finds indeed spectral statistics intermediate between Poisson and COE statistics. 

The two-point correlation function has been work out explicitly by Berkolaiko et al (2001) and has 
been shown to coincide with the statistics of so-called Seba billiards, that is, rectangular billiards with a 
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single flux line. The first few terms in a power series expansion of the form factor have been derived by 
Kottos and Smilansky (1999) and Berkolaiko and Keating (1999) and yield 



K(t) 




(15) 



see Fig. |3Ji. 



Diffusive networks: The quantum mechanics of classically diffusive systems has been studied mainly 
in the context of Anderson localisation and localisation-delocalisation transitions, see e.g. Dittrich (1996) 
and Janssen (1998) for recent review articles. 

As a simple example, we consider here a quantum graph corresponding to a classical Markov process 
on a regular lattice in d dimensions, see Fig.^t for d = 2. Choosing a stochastic matrix with constant 
transition probabilities Uj = l/2d between connected edges corresponds to rf-dimensional diffusion in the 
continuum limit L — ► oo; here, L is the number of vertices along each direction, that is, the total number 
of arcs is he = 2dL d . Solving the diffusion equation with periodic boundary conditions allows one to 
recover the low lying part of the spectrum of T^, that is, 



with diffusion constant D = ^ and m is a d-dimensional integer lattice vector. The influence of the 
dimension d on the small r behaviour of the form factor in diffusive systems has been described in detail 
by, for example, Dittrich (1996) and references therein. In terms of the spectral gap condition, one obtains 



that is, we expect universal statistics for d > 3 only. (Actually, the bound by Gnutzmann and Altland 
(2004) guarantees random matrix statistics only for d > 4, numerics suggests however, that d = 3 follows 
RMT already (Tanner 2001)). One finds convergence to the Poisson limit due to Anderson localisation in 
one dimension (Schanz and Smilansky 2000). The two-dimensional case is critical with A ~ 1/n.E which 
shows up in the form factor as a plateau for r — * 0, that is, lim r ^ -^( T ) = V^ 71 ", see Fig. |3Jd (Tanner 
2002). 

3 Regular quantum graphs 

In the following, we will consider quantum graphs for which the above bounds do not necessarily hold 
due to length correlations in the graphs, that is, averages are not taken over a full USE. We will actually 
look at a specific set of such non-generic quantum graphs, namely quantum graphs for which the global 
propagator S consists of identical local scattering processes at every vertex. Such graphs have been 
called regular quantum graphs by Severini and Tanner (2004); the name derives from the notation regular 
graphs for graphs which have the same number of incoming and outgoing arcs at every vertex and thus 
share the property that vertices are locally indistinguishable. We will show that regular quantum graphs 
corresponding to the same underlying graph G can behave very differently depending on how the local 
scattering processes are connected to each other. We will show that a crucial element in this is played by 
the possible ways regular graphs can be edge-coloured. Regular graphs are in fact another way at looking 
at quantum random walks as will be pointed out at the end of this section. 

3.1 Regular quantum graphs and edge-colouring matrices 

We will construct a quantum graph on a (i-regular digraph G with n vertices for which the wave dynamics 
at a given vertex of the graph is "locally indistinguishable" from that of any other vertex of the graph. 
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Figure 4: A 3-regular graph of size 4 together with a possible edge-colouring; the edge-colouring matrices 
correspond to entries having the same colour in the adjacency matrix A of the graph. 



This is done by choosing a unitary d-dim. matrix a and a set of arc-lengths L i: i = 1, . . . d and ascribing 
the scattering process a to every vertex in the graph with incoming as well as outgoing arcs chosen from 
the set of L;'s at every vertex. This is done here by first fixing a so-called edge- colouring of the graph 
(see eg Bollobas 1979), that is, we assign one of d different "colours" to every directed arc of the graph 
in such a way that no vertex has two incoming or two outgoing arcs of the same colour. Note that there 
are many different ways to edge-colour a given regular graph for d > 2. Edge-colouring can be described 
in terms of a set of n - dimensional permutation matrices pi,i = 1, . . .d having the property that 



E 



pi 



(17) 



with A G being the adjacency matrix of G; the arc (ij) is then assigned the colour k if the (ij) matrix 
element of pf. is non-zero. We refer to the set of pi's as the edge- colouring matrices of G . In Fig. an 
explicit example is shown with 
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The adjacency matrix A LG of the line-graph of G can then be written in the form (Severini 2003, 
Severini and Tanner 2004) 



LG 



(18) 



where 3d is the {d x d) matrix with all elements being equal to 1 and I n is the identity matrix. We 
proceed by defining a quantum graph on the line-graph of G in the form of a unitary propagator S G as 
follows: 



sG = 0^ •( C ® 1 ») = 



Cll Pi C\2 Pi 

C21 pi G22 Pn 

Cdl Pd Cd2 Pd 



CmPi 
Cid pi 

Cdd Pd 



(19) 



with C(k) = D(fc)cr and D{k)a — 8jie lkLj describing the local scattering process. The d - dimensional 
matrix C is also called the coin in the context of quantum random walks on graphs (Kempe 2003). 

Note that different ways of edge-colouring the graph, that is, different decomposition of A G in the 
form l|17J) lead to different quantum graphs which may have quite different properties as will be shown 

Severini and Tanner (2004), these matrices have been called connectivity matrices. 
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Figure 5: The complete graph K n with n = 5. 



in the next section; this is in contrast to the representations of the line-graph adjacency matrix i|18|) 
which are all equivalent up to relabeling the arcs in the graph. For modifications of this construction for 
undirected graphs with time reversal symmetry, see Severine and Tanner (2004). 



3.2 Prom 'integrable' to 'chaotic' regular quantum graphs - some examples 

In this section, we will show that different ways to edge-colour a graph can indeed lead to very different 
types of quantum graphs with spectral statistics ranging from Poisson to CUE. In other words, in regular 
quantum graphs it is the choice of the edge-colouring matrices pi, a purely topological quantity, which 
determines the properties of the quantum graph independent of the single vertex scattering processes 
given through the coin C. We will demonstrate this here for a specific example, namely so called complete 
graphs G — K n with adjacency matrix A G = J„, that is every vertex is connected to every other vertex, 
see Fig. (Note that tie — n 2 , here). A more general treatment can be found in Severini and Tanner 
(2004). In the examples discussed below, we make use of the fact that for finite groups T of order n, we 
may write 

n 

5>*=J„. (20) 

z=l 

where the p^s form a regular representation of T. In what follows we will study various decompositions 
of J„ and see how they effect statistical properties of the spectra of quantum graphs. 



3.3 The cyclic group Z n 

We will first consider an abelian group, namely the cyclic group Z„. The p^s forming a regular repre- 
sentation commute with each other and are of the form 

(Pj)ki = mod n with eigenvalues xL = e 27 " 2 ^ , where j, m = 1, . . . , n . 

The abelian nature of the group allows one to block-diagonalise the matrix S G into n blocks of dimension 
n each, independent of the coin C. The spectrum of the quantum graph is then given by the spectra of 
the sub-matrices 

S G = ^0 e 2 ^ j • C with m = 1, . . . , n . 

The eigenvalues of S G are here characterised in terms of two quantum numbers, an 'angular momentum' 
quantum number m and a second quantum number counting the eigenvalues in each m manifold. If the 
spectra for different m are uncorrclatcd, one expects Poisson statistics of the total spectrum in the limit 
n — > oo. 

Figure E^) shows spectral properties of S G with n — 24, that is, dim5 G = 576. We plot here the 
nearest neighbour spacing (NNS) distribution P(s) and the form factor K(t). The coin is of the from 



9 



C(k) = D(k) a where the local scattering matrix a is choosen to be the Fourier matrix and the arc lengths 
entering the diagonal matrix D are chosen independently and identically distributed in [0, 1]. The average 
is, for a fixed choice of the coin, taken by averaging over the wavelength k. The numerical results are 
shown in Fig. [(J,) and suggest indeed Poisson-statistics apart form deviations in the form factor on scales 
t < 1/n due to the 'random nature' of the coin. 

3.4 The non-abelian case: the symmetric group S4 

Next, we consider a specific example of a non-abelian group, namely the symmetric group S4 with n = 24 
elements. By writing the permutation matrices pi forming the regular representation of S4 in terms 
of the irreducible representations (for short irreps) of S4, we can again give the propagator S K in 
block-diagonal form where the blocks are now of size n d, with d being the dimension of the irrep under 
consideration. The group S4 has 2 one-dimensional, 1 two-dimensional and 2 three-dimensional irreps, 
and each d - dimensional irrep is contained d times in each pi according to the general formula 

2 ■ l 1 + 1 • 2 2 + 2 - 3 3 = 24. 

We thus have 5 independent sub-spectra making up the spectrum of the quantum graph S G of which two 
are of dimension n = 24, one is of dimension 2 • 24 = 48 and two are of dimension 3 • 24 = 72; the latter 
once are two and three times degenerate, respectively. The huge degeneracy in the spectra can clearly be 
seen in the spectral statistics; it is manifest in the peak at s — in P(s), see Fig. Eh, and leads to 

K(t) = (2 • 3 3 + 1 ■ 2 3 + 2 • l 3 ) = 8/3 for r>3/24. 

The spectra appear to be uncorrelated otherwise; note however, that the spectrum for each sub-block 
alone are correlated following CUE statistics, which gives rise to the deviations from purely Poisson 
behaviour in P(s) (cf. dashed curve) as well as in the behaviour of the form factor for r < 3/24 which is 
dominated by the sub-spectra of the three dimensional irreps. 

3.5 The generic case: no symmetries 

The overwhelming number of decompositions of the form H20JI will of course have no common symmetry, 
that is, it is not possible to block-diagonalise the pi's simultaneously. We therefore do not expect any 
special features in the spectrum. The question remains, however, if the 'randomness' put into the system 
by choosing a random edge-colouring is enough to produce generic, random matrix type, statistics. After 
all, these quantum graphs still possess a large degree of degeneracy due to the presence of identical coins 
at every vertex. A numerical study may thus reveal interesting insights into the range of validity of the 
RMT - regime. Fig. shows the level statistics for an unstructured choice of edge-colouring matrices 
which is in good agreement with random matrix theory for the CUE - ensemble 2 . Note that the statistic 
has been obtained from the spectrum on an (n 2 x n 2 ) unitary matrix which has only n 3 non-zero elements 
of which only n 2 are independent; in addition, there are only n different arcs lengths to choose from for 
n 2 different arcs. The origin of the universality in the spectral statistics in this type of quantum graphs 
is here clearly not due to the 'randomness' in the choice of the matrix elements but due to random 
edge-colouring of the graph alone! 

3.6 Quantum random walks 

In recent years, the study of unitary propagation on graphs has also been looked at from the perspective 
of devising a quantum version of a random walk. This line of thought arose as part of the effort to 
build quantum information systems being able to do operations which are impossible or much slower on 
classical devices. It could indeed be shown that quantum random walks have this property under certain 

2 Deviations in the formfactor for small r can be attributed to the fact that the spectrum of C itself is always contained 
in the full spectrum, see Severini and Tanner (2004); this part has been removed in the NNS statistic. 
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(a) The cyclic group (n=24): 
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1 1 1 1 1 1 1 

\ data 


0.8 
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\; Poisson 


0.6 




0.4 




0.2 




n 


1 1 1 1 1 1 
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(b) The symmetric group (n=24): 
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Poisson 
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(c) The generic case; only trivial symmetry (n=20): 
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0.5 1 1.5 2 2.5 3 3.5 4 



Figure 6: Formfactor K(t) and nearest neighbour spacing distribution P(s) for: (a) the p^s form the 
regular representation of the cyclic group Z24; (b) pi's represent the symmetric group S4; (c) a 'random' set 
of pi's without symmetries. The dashed curve in (b) labeled "red. Poisson" corresponds to a distribution 
of degenerate levels being Poisson distributed otherwise. 



11 



circumstances: quantum walks can be faster for some network geometries (Aharonov et al 1993) and can 
even lead to an exponential speed-up such as for the graph-traversal algorithm by Childs et al (2003); 
for an introductory overview and further references, see Kempe (2003). The generalisation of Grover's 
algorithm (Grover 1997) to spatial searches on graphs (Shenvi et al 2003, Ambianis et al 2004, Childs 
and Goldstone 2004) is another remarkable recent result. 

We will not study the properties of quantum random walks here; instead, we would like to point 
out that the discrete quantum walk modules discussed in the literature are in fact equivalent to regular 
quantum graphs such as introduced in the previous sections. 

In brief, a discrete quantum random walk comprises of an in general d - regular graph G with n 
vertices (where lattices are favoured in the literature) and a 'spin' degree of freedom, where the spin 
can take up d different states. A quantum state at a vertex v with spin i is then described by \v,i > 
and the quantum wave function is a superposition of all possible vertex states. The quantum random 
walk consists of a 'walk' element and a 'quantum coin toss'. The walk is steered by the internal spin 
states, that is, there exist matrices pi, i = 1, . . .d which define in which direction the component \v,i > 
is flowing, that is, 

pi\v,i >= \v',i > 

where the internal spin states remains unchanged. The pj's are obviously permutation matrices in the 
vertex space. The coin toss is simulated by applying a unitary transformation C to the spin states at 
every vertex, that is, 

C\v,i >= \v,^2,Cijj >; 
j 

the coin is unbiased, if |Cy| 2 = 1/d throughout. It is clear from the construction, that one obtains a 
classical random walk on such a network, if one performs a measurement after every walk by projecting 
out the spin degrees of freedom. 

A full cycle of walk and coin toss is then described by a unitary matrix 



qw= ©A -(C®I n ) (21) 




which has exactly the form of the propagator for regular quantum graphs in Eq. (jl9(l . We thus identify 
the p^s with the edge-colouring matrices and the coin is in fact a local scattering matrix. Writing the 
components of the wave function in \v, i > notation is thus just a different way of labeling the arcs in the 
graph and the dynamics takes indeed place on the line-digraph of G. 

To demonstrate the mechanism, let us give a specific example: the quantum walk on an infinite line. 
A popular setting is the following walk: 

states: |«,T>, \i, l>, i<El> 

walk: \i i> -> \i - 1 J> ' 1S ' ^ = i ' i+u Pi = 



n If 1 1 

com: c = 7H-i i 

A random walk starting in the state \i, |> then evolves according to 

\i, t>^ -±= T> |>) w - lk (\i + i,T> -\i - U>) ^ • • • ■ 

Continuing this process leads quickly to a probability profile which grows linearly with the number of 
steps t which is in contrast to the classical spreading for Id random walks being of the order \ft (Kempe 
2003). It must be pointed out, however, that this is a pure symmetry effect due the choice of edge- 
colouring matrices; in the setting chosen here, the two operations p\ and pi commute. By introducing 
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disorder into the quantum random walk on a line by for example randomly switching the shift operation 
at every vertex leads to a completely different behaviour resembling that of a quantum graph on a diffu- 
sive network as discussed in sec. 12.31 
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